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I. INTRODUCTION 

Gamma-ray bursts (GRBs) are the most luminous as- 
trophysical events observable today, because they are at 
cosmological distances [l| . The duration of a gamma-ray 
burst is typically a few seconds, but can range from a 
few milliseconds to several minutes. The initial burst at 
gammay-ray wavelengths is usually followed by a longer- 
lived afterglow at longer wavelengths (X-ray, ultraviolet, 
optical, infrared, and radio). Gamma-ray bursts have 
been detected by orbiting satellites about two to three 
times per week. Most observed GRBs appear to be col- 
limated emissions caused by the collapse of the core of 
a rapidly rotating, high-mass star into a black hole. A 
subclass of GRBs (the "short" bursts) appear to originate 
from a different process, the leading candidate being the 
collision of neutron stars orbiting in a binary system. See 
Ref. Q for a recent review on GRBs. 

GRBs can be used as distance indicators [s*] , thus pro- 
vide a complementary probe of dark energy^. The main 
advantage of GRBs over Type la Supernovae (SNe la) is 
that they span a much greater redshift range (from low z 
to z > 6). The main disadvantage is that GRBs have to 
be calibrated for each cosmological model tested (see for 
example, Ref. Q). This is in contrast to SNe la, where 
the calibration relations are established using nearby SNe 
la, and applied to high z SNe la to extract cosmologi- 
cal constraints. There are no nearby GRBs that can be 



used for calibration. Thus, the GRB data must be fitted 
simultaneously for calibration and cosmological parame- 
ters. This makes the use of GRBs to probe cosmology 
somewhat cumbersome. 

In this paper we show that the current GRB data can 
be summarized by a set of model-independent distance 
measurements. Our results provide an easy, robust, and 
transparent way to incorporate GRB data in an analysis 
of combined cosmological data to constrain dark energy. 



II. METHOD 

A. Calibration of GRBs 

Following we consider five calibration relations for 
GRBs. These relate GRB luminosity, L, or the total 
burst energy in the gamma rays, E^, to observables of 
the hght curves and/or spectra: riag (time lag), V (vari- 
ability), i?peak (peak of the vFn spectrum), and trt (min- 
imum rise time): 



^ Ref . _4] contains reviews with extensive lists of references on dark 
energy research. 
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Not surprisingly, i5poak carries the most distance in- 
formation. The iJpeak ~ relation is the tightest of 
the GRB calibration relations. To be included in this 
relation, the GRB afterglow must have an observed jet 
break in its light curve, and this means that only a frac- 
tion of GRBs with redshifts can contribute to establishing 
this relation. The variability-luminosity relation has the 
largest scatter. The variability is a measure of the sharp- 
ness of the pulse structure, which is determined by the 
size of the visible region in the jet. 

In order to calibrate GRBs, L and Ej must be related 
to the observed bolometric peak flux, Pboio, and the bolo- 
mctric fluence, S'boio: 



L 



47rdiPbo 



— -E7,iso^bcam — 47rd|S'bolo(l + 2) ^Fboam, (6) 

where -E-y.iso is the isotropic energy. Clearly, the calibra- 
tion of GRBs depend on the cosmological model through 
the luminosity distance dL{z). 

The cosmological constraints from GRBs are sensitive 
to how the GRBs are calibrated. Calibrating GRBs using 
Type la supernovae (SNe la) gives tighter constraints 
than calibrating GRBs internally In this paper, we 
choose to calibrate GRBs internally, without using any 
external data sets, so that our results can be used to 
combine with any other cosmological data sets. 

In fitting the five calibration relations, we need to fit 



a data array {xi,yi} with uncertainties {ax 
straight line 

y = a + bx 
through the minimization of given by 
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where we have defined 

d^={l + z)Hor{z)/c. 



(15) 



Since the absolute calibration of the GRBs is unknown, 
the Hubble constant cannot be derived from GRB data. 
Thus we have defined the data arrays {yi} such that 



c/Ho = 9.2503 X lO^"^/!"^ cm 
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21ogdL, 



(14) 



is absorbed into the overall calibration. 

Furthermore, for the L - E-p^ak and i?pcak - E^ re- 
lations, the measurement error of iJpoak is asymmetric, 
thus we need to modify the such that 

= cr+j, if (2/i - a)/6 > a;^; 

Cfx.i = if {yi - a)/b < Xi, (17) 

where cr^,- and a~ ^ are the ± measurement errors. 

As noted by Ref.Q, the statistical errors on {ai,bi} 
are quite small, but the x^^s are very large due to the 
domination of systematic errors. Following Ref.0], we 
derive the systematic errors by requiring that — 
(the degrees of freedom), and that a^^^ = a'^f^t + '^'sys- 

For illustration on the cosmological parameter depen- 
dence of the calibration of GRBs, Table 1 shows the sys- 
tematic errors, as well as the constants {ai,bi} for the 
five calibration relations for =0.2, 0.27, 0.4 for a flat 
Universe with a cosmological constant. Note that the 
in this table is smaller than the definition of Ref.[6j] by 
21og(9.2503 X lO'^'^h-^). Note also that the derived sys- 
tematic errors change by less than 3% for the different 
models. Since systematic errors should be independent 
of the cosmological model, we take the systematic errors 
for the rim = 0.27 flat ACDM model to be the standard 
values in the rest of this paper. 

For most of the calibration relations, fitting straight 
lines using errors in both coordinates does not give sig- 
nificantly different results from fitting straight lines using 
errors in the y coordinates only. For the L - V rela- 
tion, assuming the flm = 0.27 fiat ACDM model, fit- 
ting straight lines using errors in the y coordinates only, 
we find 02 = -3.540 ± 0.003, 62 = 1-649 ± 0.012, and 
<^sys,2 = 0.518. The slope and the systematic error are 
both significantly smaller than the results shown in Table 
1. Since the x coordinates have significant measurement 
errors in all five calibration relations, the latter should 
be fitted to straight lines using errors in both x and y 
coordinates. 

Note that we do not use the calibration relations from 
Table 1 (the parameters and bi) when we derive model- 
independent distances from GRBs (see Sec.IIB). Table 1 
is only used to show that the calibration of GRBs is sensi- 
tive to the assumptions about cosmological parameters, 
but the systematic uncertainties of the calibration pa- 
rameters ai and bi are not sensitive to the assumptions 
about cosmological parameters. Thus, we will derive cal- 
ibration relations of GRBs for each set of assumed dis- 
tances, but we will assume that the systematic errors on 
tti and bi are given by the = 0.27 flat ACDM model. 
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TABLE I: Systematic errors for the five GRB calibration relations 
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B. Model-independent distance Measurements 
from GRBs 

Following Ref.fgl, we weight the five estimators of dis- 
tance of GRBs (from the five calibration relations) as 
follows 



(log^i). ^ — A. . /: — - (18) 



, data 



2\ data 



-1/2 



(19) 



where 



2\ (^) ( \ 



for a — 1,2,3,5 



(20) 



/ 2\(4) (4) 

{log dh j = 04 + bix] - log 



47rS'bolo,i^be 



and 



l + z 



(21) 



,4"^nlO, 



0'(-Pbolo,i) 



.^'boio.ilnlO 
for a — 1, 2, 3, 5 



\ sys j 



(22) 



Ultimo. 



bolo.-i ) 



•Sboio.i In 10 



*^(-^bcam,i) 



^bc 



I In 10 . 



(23) 



For each cosmological model given by d^ , we can cali- 
brate the GRBs as described in Sec lII Al and then derive 

the distance estimate (log di,^^ from each GRB. The 

of a model is given by 
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The treatment of the asymmetric errors in -Bpcak is given 
by Eq.ini). 

We summarize the cosmological constraints from GRB 

data, (zgrb), in terms of a set of model-independent 
distance measurements {r^{zi)}: 



rp{zi) ^ 



<i), 

rp(0.l7) z cn 



(25) 

where r{z) = dL{z)/{l + z) is the comoving distance at z. 
For the 69 GRBs from [Q], the lowest redshift GRB has 
z — 0.17, while the highest redshift GRB has z = 6.6. 
There are only four GRBs at 4.5 < z < 6.6. We find 
that the optimal binning is to divide the redshift range 
between 0.17 and 4.5 into 5 bins, and choosing the last 
bin to span from 4.5 and 6.6 (see Fig. 2). 

Note that the ratio rp{z)/rp{0. 17) is the most conve- 
nient distance parameter choice for the currently avail- 
able GRB data, since z = 0.17 is the lowest redshift GRB 
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FIG. 1: The accuracy with which the scaled distance T^iz) = 
rp(2:)/rp(0.17) can be reconstructed using cubic spline interpolation 
(C.S.) from {?p(^i)} for Afj,i„ = 6 with our choice of binning. Note 
that Af7(2) = T^{z)c.s. -r^{z)true- 



in the data set, and the absolute cahbration of GRBs 
is unknown. Using the distance ratio rp{z)/rp{0.17) 
removes the dependence on Hubble constant (which is 
unknown due to the unknown absolute calibration of 
GRBs). 

Because {r^{zi)} varies very slowly for all cosmological 
models allowed by current data, the scaled distance r^{z) 
at an arbitary redshift z can be found using cubic spline 
interpolation from {r^(zi)} to ~ 1-3% percent accuracy 
for Nbin = 6 with our choice of binning, see Fig[TJ 

Note that for a given set of possible values of {r^{zi)} 
[i = 1,2, ...,6), the luminosity distance at an arbitrary 
redshift, d^^z), is given by the accurate interpolation de- 
scribed above. Thus no assumptions about cosmological 
parameters are made. We calibrate the GRBs for each set 
of possible values of {r^{zi)} (i = 1, 2, 6), and compute 
the likelihood of this set of {r^{zi)} in a Markov Chain 
Monte Carlo analysis. Hence the distances {r^{zi)} are 
independent of assumptions about cosmological parame- 
ters. 



C. Other cosmological data 



GRB data, we study them in combination with 307 SNe 
la Cosmic Microwave Background anisotropy (CMB) 
data from WMAP five year observations [l3| , and baryon 
acoustic oscillation (BAO) scale measurement from the 
SDSS data [llj. 

SN la data give the luminosity distance as a function 
of redshift, dL{z) = {1 + z) r{z). We use 307 SNe la from 
the "Union" compilation by Ref.^[jWhich includes data 
from Ref. [l2| • Appendix A of Ref. l3| describes in detail 
how we use SN la data (flux-averaged to reduce lensing- 
like systematic effects [1^ and marginalized over Hq) in 
this paper. ^ We applied flux-averaging to the "without 
systematics" data from the "Union" compilation. The 
"with systematics" data differ only in having small off- 
sets added to different data sets, which leads to correla- 
tions in the data. We expect the increase of uncertainties 
resulting from flux- averaging to be larger than the effect 
of the small offsets between different data sets. 

We include the CMB data using the method proposed 
by Ref. 14], which showed that the CMB shift parameters 
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(26) 



together with fli,h^, provide an efficient summary of 
CMB data as far as dark energy constraints go. We 
use the covariance matrix of [R{z^),la(z^:), rihh^] from the 
five year WMAP data (Table 11 of 0), with z^given by 
fitting formulae from Hu & Sugiyama (1996) [i5|. CMB 
data are included in our analysis by adding the follow- 
ing term to the of a given model with pi — R{z^:), 
P2 = laiz*), and P3 = flbh'^: 

XcMB = [Cov^\jg{p^,pj)] Apj, Ap, = Pi 



•Pi 



^data 
i 5 

(27) 

where are the maximum likelyhood values given in 
Table 10 of [13. 

We also use the SDSS baryon acoustic oscillation 
(BAO) scale measurement by adding the following term 
to the of EL model: 



Xbao 
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where A is defined as 
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(29) 



and Abao = 0.469 (ns/0.98)-*'-35, a a = 0.017, and 
zbao = 0.35 (independent of a dark energy model) [ll|. 
We take the scalar spectral index rig = 0.96 as measured 
by WMAP five year observations ^101] . 

Finally, in combination with CMB data, we include 
the Hubble Space Telescope (HST) prior on the Hubble 
constant of = 0.72 ± 0.08 [16|]. 



GRB data alone do not constrain dark energy pa- 
rameters. In order to investigate how well our model- ^ A public code for fl ux-averaging SN la data is available at 

independent distance measurements from GRB represent |http://www.nhn.ou!edu/~wang/SNcode/| 
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III. RESULTS 



For Gaussian distributed measurements, the likelihood 
function L oc e~-^ ' , with 

= Xgrb + xInc + XcMB + Xbao^ (30) 

where Xgrb given in Eqs. ((3T)) - ((33)l . XsNe given in 
Appendix A of Ref. , Xcmb i^ given in Eq. ([27|) , and 
xIao is given in F.q.^. 

We run a Monte Carlo Markov Chain (MCMC) based 
on the MCMC engine of [13 to obtain O(IO^) samples 
for each set of results presented in this paper. For the 
model-independent GRB distance measurements, the pa- 
rameter set is {r^{zi)} {i = 1, 2, 6); no assumptions are 
made about cosmological parameters. For the combined 
analysis of GRBs (using either the model-independent 
distance measurements {r^{zi)} {i — 1, 2, 6), or the 69 
GRBs directly) with other cosmological data, the cosmo- 
logical parameter sets used are: for a flat Universe 
with a cosmological constant; (rim, ^^a) for a cosmologi- 
cal constant; (Jim, wq) for a flat Universe with a constant 
dark eneagy equation of state; and (fim, h, Vl},h?, Pde) 
for a flat Universe with a dark eneagy equation of state 
linear in the cosmic scale factor, with pde = (wo, W0.5) 
or (uiQ, Wa)- We assum flat priors for all the parameters, 
and allow ranges of the parameters wide enough such 
that further increasing the allowed ranges has no impact 
on the results. The chains typically have worst e-values 
(the variance(mean)/mean(variance) of 1/2 chains) much 
smaller than 0.005, indicating convergence. The chains 
are subsequently appropriately thinned to ensure inde- 
pendent samples. 

FigO shows the distances {r^{zi)} measured from 69 
GRBs using the five calibration relations in Eqs.(IT])-([5]). 
Table 2 gives the mean and 68% confidence level errors of 
\jp{zi)}- The normalized covariance matrix of {r^{zi)} is 
given in Table 3. To use our GRB distance measurements 
to constrain cosmological models, use 

Xgrb = [Ar>(zi)] ' (Covq^^) • [A7>(zj)] 
Arv(z,) = r>'^"*"(z,)-7>(z.), (31) 

where 7^(2) is defined by Eq. (|25p . The covariance matrix 
is given by 

(CovGfls)„- = cr(7>(z,))cr(r>(2;j)) (CovGits)^^. , (32) 

where Covgrb is the normalized covariance matrix from 
Table 3, and 

(T(rv(z,)) - a (r^(z,))+ , if r^iz) > 7>(z)'^^*^ 
a(rv(z.)) - a(7v(z.))-, if 7^ W < 7v(z)^'^'%(33) 

where a (r^{zi))^ and a (j^^Zi))" are the 68% C.L. errors 
given in Table 2. 

Using the distance measurements from GRBs (see Ta- 
bles 2 and 3), we find fl^ =0.247 (0.122, 0.372) (mean 
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FIG. 2: The distances measured from 69 GRBs using the five 
calibration relations in Eqs.JTJ-©. The error bars indicate the 
68% C.L. uncertainties. The top panel shows the scaled distances 
rp{zi) (see Ea. ll25l l): the bottom panel shows the corresponding 
comoving distances r(zi). 



and 68% C.L. range). Using the 69 GRBs directly, 
we find =0.251 (0.135, 0.365). This demonstrates 
that our model-independent distance measurements from 
GRBs can be used as a useful summary of the current 
GRB data. 

FigOshows the joint confidence contours for ($!,„ , JIa), 
from an analysis of 307 SNe la with and without 69 
GRBs, assuming a cosmological constant. This shows 
that the addition of GRB data significantly reduces the 
uncertainties in {flm, ^a), and shifts the bestfit parame- 
ter values towards a lower matter density Universe. Fig[4] 
shows the joint confidence contours for {wq, fim), from 
an analysis of 307 SNe la with and without 69 GRBs, 
assuming a flat Universe. This shows that SNe la alone 
rules out a cosmological constant at greater than 68% 
C.L., but the addition of GRB data significantly shifts 
the bestfit parameter values, and a cosmological constant 
is consistent with combined SN la and GRB data at 68% 
C.L. 

We now consider a dark energy equation of state linear 
in the cosmic scale factor a. FigO shows the joint confi- 
dence contours for (wq, W0.5) and (wq, Wa), from a joint 
analysis of 307 SNe la with CMB data from WMAP5, 
and SDSS BAO scale measurement, with and without 69 
GRBs. HST prior on Hq has been imposed and a flat 
Universe is assumed. Note that wq.s = wx{z = 0.5) in 
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TABLE II: Distances measured from 69 GRBs with 68% C.L. upper and lower uncertainties. 
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TABLE III: Normalized covariance matrix of distances measured from 69 GRBs 
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FIG. 3: The joint confidence contours for {flm, ^a), from an 
analysis of 307 SNe la with and without 69 GRBs. A cosmological 
constant is assumed. 



the linear parametrization [16 
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with flc — 2/3 (i.e., Zc = 0.5). Eq. (|M)l corresponds to a 



dark energy density function 
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exp < 3 
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ln(l 
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Eq. ([M|) is related to wx{z) 
setting [18] 



Wo + {I - a)wa 



(35) 
by 



wo.5 - Wq 

1 - flc 



or Wo. 5 = Wo + (1 - ac)Wa. 

(36) 

Ref.[18l| showed that {wq, W0.5) are much less correlated 
than (wO, Wa), thus are a better set of parameters to 
use. FigO shows that the addition of GRB data notably 
shifts the 68% C.L. contours of (wqjWq.s) and {wq, Wa) to 
enclose the cosmological constant model {wx{a) = — !)• 

IV. SUMMARY AND DISCUSSION 

We have shown that the current GRB data, consist- 
ing of 69 GRBs spanning the redshift range from 0.17 
to 6.6, can be summarized by a set of distance measure- 
ments (see Figl2]). For each set of possible distance val- 
ues, the GRBs are calibrated. The resultant distance 
measurements (given in Tables 2 and 3) are independent 
of cosmology, and can be easily used to combine with 
other cosmological data sets to constrain dark energy (see 
Eqs.[SI]-IS3] and [15]). 

The number of bins used in our distance measurement, 
Tibin = 6, is determined by the current sample of GRBs. 
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FIG. 4: The joint confidence contours for (wq, Qm), from an anal- 
ysis of 307 SNe la with and without 69 GRBs. A fiat Universe and 
dark energy with constant equation of state are assumed. 
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Increasing the number of bins by one leads to oscillations 
in the measured distance values. As more GRB data 
become available, we can expect to be able to increase 
the number of bins used to represent GRB data. 

We find that GRB data alone give f2„ = 0.25lo;ii for 
a flat Universe with a cosmological constant. Figl3]and 
FigH] show that combining GRB data with SN la data 
significantly shifts the bestfit model towards a lower mat- 
ter density Universe, in agreement with galaxy redshift 
survey [lol and CMB data [2l[. Fig [5] shows that as- 
suming a dark energy equation of state linear in cosmic 
scale factor a, including GRB data together with SN la, 
CMB, and BAG data shifts the 68% G.L. contours of the 
two dark energy parameters to enclose the cosmological 
constant model. 

For a flat Universe with a cosmological constant, our 
results for GRBs alone differ from that of Ref.Q. We 
find f7,„ = 0.25to \l, while Ref.jl] found = 0.39+!^;^^. 
The difference likely resulted from the different numer- 
ical methods used to calibrate the GRBs. We have fit- 
ted straight lines to the calibration relations using errors 
in both X and y coordinates (see Eq.[S]), while Ref. 
used "ordinary least squares without any weighting'' . 
Our results are consistent with the results of Ref.[3], 
who used SNe la to help calibrate the GRBs and found 

^^ni — — 0.05' 

In our quest to solve the mystery of dark energy, GRBs 
will provide a unique and complementary probe. Our 
results will make it very convenient to incorporate GRB 
data in any joint cosmological data analysis in a simple 
and robust manner. 



FIG. 5: The joint confidence contours for {wo, wo.a) and {wo, Wa), 
from a joint analysis of 307 SNe la with CMB data from WMAP5, 
and SDSS BAO scale measurement, with and without 69 GRBs. 
HST prior on Hq has been imposed and a flat Universe is assumed. 



Acknowledgements I acknowledge the use of getdist 
from cosmomc in processing the MCMC chains. 



[1] Paczynski, B., 1995, PASP, 107, 1167. 

[2] Meszaros, P., Rep. Prog. Phys., 69, 2259 (2006) 

[3] Ghirlanda, G. et al. 2004, ApJ, 613, LIS; Dai, Z. G., 
Liang, E. W., & Xu, D. 2004, ApJ, 612, LlOl; Fir- 
mani, C., Ghisellini, G., Ghirlanda, G., & Avila-Reese, 
V. 2005, MNRAS, 360, LI; Liang, E. W. & Zhang, B. 
2005, ApJ, 633, 603; Firmani, C., Avila-Reese, V., Ghis- 
ellini, G., & Ghirlanda, G. 2006, MNRAS, 372, L28; 
Ghirlanda, G., Ghisellini, G., & Firmani, C. 2006, New 
J. Phys, 8, 123; Firmani, C., Avila-Reese, V., Ghisellini, 
G., & Ghirlanda, G. 2007, RMxAA, 43, 203; Kodama, 
Y. et al. 2008, MNRAS in press (arXiv:0802.3428 ); Am- 
ati, L. et al. 2008, arXiv:0805.0377; Basilakos, S., & 
Perivolaropoulos, L. 2008, arXiv:0805.0875 , Capozziello, 
S. & Izzo, L. |a rXiv:0806.1120 Wei, H., & Zhang, S. N. 
2008, arXiv:080 8.2240 Liang, N., & Zhang, S. N. 2008, 
irXiv:0808.2655] 

[4] Padmanabhan, T, Phys.Rep., 380, 235 (2003); Peebles, 
P. J. E., & Ratra, B., Rev. Mod. Phys., 75, 55 (2003); 



Sahni, V., & Starobinsky, A., IJMPD, 15, 2105 (2006) 
Gopeland, E. J., Sami, M., & Tsujikawa, S., IJMPD, 15, 
1753 (2006); Ruiz-Lapuente, P., Class. Quantum. Grav., 
24, 91 (2007); Ratra, B., & Vogeley, M. S., PASP, 120, 
235 (2008) 

[5] Schaefer, B. E. 2004, ApJ, 602, 306. 

[6] Schaefer, B. E., 2007, ApJ, 660, 16 

[7] Liang, N., Xiao, W. K., Liu, Y., Zhang, S. N. 2008, ApJ, 
in press, arXiv:0802.4262 

[8] Press, W.H., Teukolsky, S.A., Vettering, W.T., & Flan- 
nery, B.P. 2007, Numerical Recipes 3rd Edition, Cam- 
bridge University Press, Cambridge. 

[9] Kowalski, M., et al., 2008, arXiv;0804.4142 
[10] Komatsu, E., et al., arXiv;0803.0547 [astro-ph] 
[11] Eisenstein, D., et al., ApJ, 633, 560 
[12] M. Hamuy, et al., AJ, 112, 2398 (1996); Riess, A. G., 
et al., AJ, 116, 1009 (1998); Riess, A. G., et al., AJ, 
117, 707 (1999); Perlmutter, S., et aL, ApJ, 517, 565 
(1999); Tonry, J. L., et al., ApJ, 594, 1 (2003); Knop, 



8 



R. A., ApJ, 598, 102 (2003); Krisciunas, K., et al., AJ, 
127, 1664 (2004a); Krisciunas, K., et al., AJ, 128, 3034 
(2004b); Barris, B. J., ApJ, 602, 571 (2004); Jha, S., et 
al., AJ, 131, 527 (2006); Astier, P., A&A, 447, 31 (2006); 
Riess, A. G., et al., ApJ, 659, 98 (2007) Miknaitis, G., et 
al., ApJ, 666, 674 (2007) 
[13] Wang, Y., ApJ 536, 531 (2000b); Wang, Y., & Mukher- 
jee, P. 2004, ApJ, 606, 654; Wang, Y., JCAP, 03, 005 
(2005) 

[14] Wang, Y., & Mukherjee, P., PRD, 76, 103533 (2007) 



[15] Hu, W., & Sugiyama, N. 1996, ApJ, 471, 542 
[16] Freedman, W. L., et al. 2001, ApJ, 553, 47 
[17] Lewis, A., & Bridle, S. 2002, PRD, 66, 103511 
[18] Wang, Y., PRD, 77, 123525 (2008) 

[19] Chevallier, M., & Polarski, D. 2001, Int. J. Mod. Phys. 
DIO, 213 

[20] Tegmark, M., et al. 2006 , Phys. Rev. D 74, 123507 
[21] Dunkley, J., et al. 2008, larXiv:0803T0586l 



